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Abstract 

We report detailed atomistic simulations combined with high-fidelity conductance calcu- 
lations to probe the structural origins of conductance fluctuations in thermally evolving Au- 
benzene- 1 ,4-dithiolate-Au junctions. We compare the behavior of structurally ideal junctions 
(electrodes with flat surfaces) to structurally realistic, experimentally representative junctions 
resulting from break junction simulations. The enhanced mobility of metal atoms in struc- 
turally realistic junctions results in significant changes to the magnitude and origin of the 
conductance fluctuations. Fluctuations are larger by a factor of 2-3 in realistic junctions com- 
pared to ideal junctions. Moreover, injunctions with highly deformed electrodes, the conduc- 
tance fluctuations arise primarily from changes in the Au geometry, in contrast to results for 
junctions with non-deformed electrodes, where the conductance fluctuations are dominated by 
changes in the molecule geometry. These results provide important guidance to experimental- 
ists developing strategies to control molecular conductance for device applications, and also to 
theoreticians invoking simplified structural models of junctions to predict their behavior. 
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Conductance fluctuations in molecular transport junctions are a major barrier to the construc- 
tion of reliable molecular-based circuitry. 1 2 The fluctuations arise from changes in the junction 
structure between successive junction rupture and reformation events or due to thermal motion.^ 
Therefore, developing strategies to suppress conductance fluctuations relies critically on under- 
standing the structural origins of the fluctuations. For instance, following a study^that showed 
that the conductance through biphenyl molecular wires depends on the dihedral angle between 
the phenyl rings, Kiguchi et al.^ synthesized a rotaxane structure to limit changes in the dihedral 
angle, thereby suppressing the conductance fluctuations. For simpler molecules, such as benzene- 
1,4-dithiolate (BDT), the conductance fluctuations are typically attributed to changes in the metal- 
molecule contact geometry (bonding site and tilt angle) J*"^ However, the electrode geometry may 
play an increasingly important role for systems involving mechanical elongation and deformation 
of the junctionJ^HSMmpOT example, Au-thiolate bonding results in significant deformation of the 
electrodes in break junction experiments. Several groups^^l have recently investigated the role 
of Au-thiolate bonding in break junction environments, but the exact structural origins of the con- 
ductance behavior remains unclear. Many computational studies have been performed to elucidate 
these features,^^ however, almost all have considered idealized electode geometries which may 
not be fully representative of the experimental counterparts. Here, we compare Au-BDT-Au junc- 
tions with structurually ideal vs. realistic features to understand (1) the origins of conductance 
fluctuations and (2) how ideality may influence the conclusions drawn from computation. 

To evaluate the effect of junction structure on conductance, we perform atomistic simulations 
combined with conductance calculations; in this approach, snapshots are periodically extracted 
from the simulations and then used as input in electron transport calculations. Previous studies 11 ^ 1 
of this kind focused on ideal geometries with a single molecule sandwiched between two flat 
surfaces. Here, we focus on geometries that are more representative of those likely to appear in 
widely used break junction experiments JIHEEHlillMIS] j n addition to simulations where all atoms 
are dynamic, we perform two separate simulations (from the same starting point) for each junction 
where either the BDT geometry or Au geometry is fixed. The BDT geometry (i.e., intramolecular 
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geometry and Au-BDT contact geometry) is fixed by treating the molecule (including the Au atoms 
covalently linked to the BDT) as a rigid body. In separate simulations, the positions of the Au atoms 
are fixed while the BDT is free to move. By eliminating specific degrees of freedom within our 
simulations, we are able to determine the independent contributions of changes in the Au and BDT 
geometries to the conductance fluctuations. 
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Figure 1: Calculated conductance histograms. (Top row) Ideal, flat-surface junction, (middle row) 
a junction with curved tips, and (bottom row) a highly deformed junction. For each junction, 
three separate simulations are run: (left column of plots) one where all atoms in the junction are 
dynamic, (middle column of plots) one with the Au atomic positions fixed, and (right column of 
plots) one with fixed BDT geometry. The standard deviation, a, is shown with each histogram, 
and the RMSD of the Au atom bonded with BDT in each tip is shown on the far right. 



Figure 1 compares the conductance histograms for the three junctions and three simulation 
types. We first consider the fully dynamic simulation results, which reveal important differences 
between the ideal and non-ideal junctions. The ideal junction (junction 1) produces conductance 
histograms that are much narrower than those for the non-ideal junctions (junctions 2 and 3). The 
peak width (standard deviation, a) is more than an order of magnitude lower than the average 
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conductance for the ideal junction, while for the non-ideal junctions the peak width is on the same 
order of magnitude as the conductance values themselves. The large peak widths in the non-ideal 
junctions result from increased geometric freedom, and may present challenges for applications 
where a device is required to maintain a target conductance value within some threshold. The shape 
of the conductance histogram also changes for the non-ideal junctions. While the distribution for 
the ideal junction appears Gaussian (as expected for nonresonant tunneling through molecules^), 
results for the non-ideal junctions exhibit long tails spanning conductance values much higher 
than the peak values. This results from a transition in the electron transport mechanism from far 
off resonance in ideal junctions to off resonance in non-ideal junctions (see Supporting Information 
for a comparison of the transmission curves). 

It is apparent from the relative peak widths in the fixed BDT results that the role of the elec- 
trodes becomes increasingly important as the electrodes are deformed. For the ideal junction, the 
Au atoms are closely bound to their lattice sites and thus do not contribute significantly to the 
conductance fluctuations. In fact, the fully dynamic peak width is completely resolved in the fixed 
Au simulation, indicating that the conductance fluctuations are dominated by the ability of BDT to 
explore configuration space. The situation changes for junction 2, as the fully dynamic peak width 
is not completely resolved from the fixed Au simulation. This suggests that the motion of the elec- 
trodes facilitates the sampling of a greater range of contact geometries. For junction 3, the peak 
widths are similar for the three types of simulations, albeit slightly wider for the fully dynamic 
simulation, which demonstrates the importance of the interplay between Au and BDT geometry in 
these systems; that is, changes in the BDT geometry are often enabled by changes in the electrode 
geometry, and vice versa. Importantly, the peak width in the fixed BDT simulation is larger than 
that for the fixed Au simulation, indicating a transition in the primary origin of conductance fluc- 
tuations from changes in the molecule geometry to changes in the Au geometry. We attribute this 
transition to the enhanced dynamic structural fluxionality 21 (lengthening and weakening of Au-Au 
bonds) in the top tip of junction 3. To support this explanation, we calculate the root-mean-square 
deviation (RMSD) of the position of the Au atom bonded to BDT relative to its average position. 
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In cases where multiple Au atoms are bonded to BDT, the Au atom that is on average closest to 
the bonded S atom is considered. As shown in Figure 1, the peak widths scale with the RMSD 
magnitude. With the BDT geometry fixed, the top tip injunction 3 exhibits the highest RMSD due 
to its low coordination. In the fully dynamic simulations, junction 2 produces the highest RMSD, 
which is consistent with its large spread in conductance. The high RMSD results from electrode 
rearrangements, contributing to the mobility beyond simple fluctuations about a single position. 
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Figure 2: Tilt angle (°) distribution during the fixed Au simulations. 



The sampled molecular tilt angles also change between the different junction geometries. How- 
ever, increases in the range of sampled tilt angles do not increase conductance fluctuations signifi- 
cantly. Figure 2 plots histograms of the tilt angle (angle between the S-S vector and z-axis) during 
the fixed Au simulations. Note that the distribution is much wider injunction 3 where the molecule 
can more easily rotate around the "sharp" upper tip. It has been shown that in junctions where the 
electrodes are represented as flat surfaces, the conductance is sensitive to the tilt angle at values 
greater than 20°.^ It is therefore surprising that the increased tilting freedom of the BDT molecule 
injunction 3 results in conductance fluctuations that are slightly smaller than those resulting from 
the fixed BDT simulation. In the case of sharp tips, the strong relationship between tilt angle and 
conductance may not apply since the interactions between the molecule and electrode(s) are lim- 



ited by the small number of metal atoms in the vicinity of the metal-molecule bond(s). In different 
environments (e.g., higher tilt angles and relatively flat tips) where the carbon atoms in BDT can 
interact with the Au tips, the range in sampled tilt angles may make more significant contributions 
to the conductance fluctuations. 
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Figure 3: (Top) Thermally averaged and single-geometry conductance trace for Au-BDT-Au junc- 
tion undergoing elongation. The initial junction geometries are shown above for every A of elon- 
gation. (Bottom) Plot showing the average RMSD of the Au atom bonded to BDT. 

An important consideration in molecular break junction experiments is how the conductance 
fluctuations change as a junction is elongated. Thus, we next explore conductance fluctuations in a 
junction undergoing mechanical elongation. The junction elongation procedure is described in the 
Computational Methods section. The average conductance and standard deviation at each elonga- 
tion length are presented in Figure 3, along with the RMSD value of the BDT-bonded Au atom in 
each tip during the MD-MC simulation; the conductance from the initial geometry at each elon- 
gation length is also plotted for comparison. The error bar sizes in Figure 3 strongly indicate that 
mechanical deformation of the junction increases the conductance fluctuations, as the error bars at 
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the early stages 0.0 and 0.5 A are very small, before increasing significantly at 1 A. This behav- 
ior is dictated by high-mobility structures that form in response to mechanical elongation; that is, 
the tips become less ideal as the junction is deformed. The fluctuations and RMSD are largest at 
elongation lengths of 1.0 and 2.0 A, where significant structural rearrangements of the electrodes 
occur during the simulation. Figure 3 also highlights the importance of considering more than a 
single geometry when calculating conductance, especially for systems where the relative changes 
in conductance are small between the various structures. It is clear that the single-geometry data 
fails to capture the average behavior over the entire range. Additionally, Figure 3 illustrates the 
difficulty in identifying junction structure based on experimental conductance histograms, as dif- 
ferent junction structures may have similar average conductance values and highly overlapping 
distributions. 

In summary, we have investigated conductance fluctuations in structurally distinct Au-BDT-Au 
junctions. We demonstrated that conductance fluctuations in non-ideal junctions are higher than 
those in an ideal junction. We also showed that while changes in the molecule geometry dominate 
conductance fluctuations in structurally ideal junctions, the enhanced motion of the Au atoms 
in deformed electrodes leads to an increased contribution to the conductance fluctuations from 
changes in the electrode geometry. The minimal role of Au geometry in flat-surface junctions and 
its significant role in highly deformed junctions highlight the importance of controlling structure in 
single-molecule conductance measurements. These results also show that conductance fluctuations 
in thiolate-based break junctions, where significant deformation to the electrodes occurs, ^ 1 ^ 13 l 
may be difficult to control as both the molecule and electrode motion make significant contributions 
to the fluctuations. The presence of other complex bonding arrangements (e.g., Au-S-Au-S-Au 
"staple" motifs^ at the Au-S interface may further complicate this issue.^^Hln contrast, linkers 
with weaker coupling (e.g., amines) are unlikely to result in the formation of structures such as Au- 
Au2-Au, and thus efforts to control molecular motion 5 may prove highly effective for controlling 
conductance fluctuations. This work also demonstrates the importance of the choice of electrode 
used in computational studies. 
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Computational Methods 

Hybrid Molecular Dynamics-Monte Carlo Simulations. In this study, we simulate three struc- 
turally distinct junctions: BDT connected between two perfectly flat Au(100) surfaces (junction 1), 
two curved tips (junction 2), and a curved tip and a Au-Au2-Au structure (junction 3). Junctions 
2 and 3 were constructed via a hybrid molecular dynamics-Monte Carlo (MD/MC) simulation ap- 
proachJ^ESI an( j were selected as they represent two extremes of the types of junctions that would 
be expected in real break junction experiments. Note, the stability and appearance of the Au-Au2- 
Au structure and other similar low-coordination structures are supported by density functional 
theory calculations J 2 ^ 2 ^ Each junction is evolved at 77 K (a common temperature in experimentP 
and simulation^ by performing 200 cycles of MD-MC sampling, where a cycle consists of 0.2 
ns of MD followed by 200,000 MC moves; at the end of each cycle, the conductance is com- 
puted. The hybrid MD-MC protocol yields ergodic sampling of the possible configurations in each 
junction; that is, the conductance histograms are converged, implying that we have fully sampled 
configuration space. In the junction elongation simulations (see Figure 3), elongation is carried out 
by displacing the upper lead by 0.1 A in the [001] direction followed by 20 ps of MD and 100,000 
MC moves. Every 0.5 A of elongation a snapshot is extracted and evolved without stretching for 20 
MD-MC cycles, with the conductance computed after each cycle. MD simulations are performed 
in LAMMPS 29 while the MC moves, which are included to enhance the sampling of the preferred 
Au-S bonding geometries, are performed using an in-house code. The second-moment approxima- 
tion to the tight-binding potential (TB-SMA) is used to describe Au-Au interactions; TB-SMA is 
a semiempirical, many-body potential that was developed by Cleri and Rosato 30 to capture metal- 
lic bonding effects. Multi-site Morse bonding curves calibrated from density functional theory 
calculations are employed to describe S-Au bonding. 31 Van der Waals interactions between BDT 
and Au and intramolecular interactions within BDT are described using the universal force field.^ 
Electrostatic interactions are described using the Coulombic potential, with partial charges resid- 
ing on BDT atoms derived from our prior work. 33 Further details are included in the Supporting 
Information. 
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Conductance Calculations. Conductance calculations are performed using the geometries ex- 
tracted directly from the hybrid MD/MC simulations; no geometry optimizations are performed 
within the DFT framework. The zero-bias conductance, given by G = T(Ef)Gq, where T is the 
transmission probability, Ep is the Fermi level of the electrodes, and Go = 2e 2 //z, is calculated us- 
ing SMEAGOLjSESIan electron transport code that interfaces with the DFT package SIESTA.E^EZl 
Self-interaction corrected DFT is employed to accurately describe the energy level lineup between 
the molecule and leads .1121211 We obtain a conductance value of O.O6G0 after optimizing the geome- 
try of junction 1 (using the classical force fields), which matches the value reported in prior studies 
of ideal Au-BDT-Au junctions,! 19 l 38 * thus validating our force fields and the conductance calcula- 
tions. This value of conductance (O.O6G0) along with most of the other values shown in Figure 1 
are higher than the most-probable value (0.01 IGo) reported in some experiments J2~E3 while falling 
in the range of values measured in other experiments.^^ As pointed out in several previous stud- 
ies J22KD the long-standing discrepancy between the experimentally measured and theoretically 
calculated value may be due to shortcomings in the theory, structural differences between the 
experimental and theoretical junctions, or some other factor. Nevertheless, these differences are 
immaterial for the present work since the focus is on relative conductance differences rather than 
the exact values of conductance. Further details are included in the Supporting Information. 

Note that we do not consider vibrational effects in the conductance calculations, instead focus- 
ing only on elastic transport processes. It has been shown 42 that vibrational effects are negligible 
for BDT connected to Au tips, and for Au nanowires vibrational effects are small compared to the 
total transmission. 43 We therefore expect a small vibrational contribution to the transport for our 
Au-BDT-Au structures. 
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